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A new kind of delta expansion is applied on the lattice to the d = 2 non-linear a model at N = oo 
and N — 1 which corresponds to the Ising model. We introduce the parameter 8 for the dilation of 
the scaling region of the model with the replacement of the lattice spacing a to (1 — <5) 1//2 a. Then, we 
demonstrate that the expansion in S admits an approximation of the scaling behavior of the model 
at both limits of N from the information at a large lattice spacing a. 
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I. INTRODUCTION 

It is one of the outstanding problems to construct a systematic computational framework to study nonpcrturbative 
aspects of quantum fields. Lattice field theories initiated by Wilson [lj allow us to use the strong coupling expansion 
which is mathematically equivalent with the high temperature expansion in condensed matter physics. However, the 
results do not necessarily provide us, at least in the quantitative sense, the corresponding results in continuum space- 
time since the strong coupling expansion on the lattice usually breaks down at small lattice spacings. Nevertheless 
the fact that the strong coupling expansion on the lattice clarified various nonperturbative properties of quantum 
fields such as the quark confinement serves us enough motivation to investigate the possibilities of improving it as to 
be effective on the approximation of physics in the continuum limit. 

As an attempt toward the improvement of the strong coupling expansion on the lattice, we propose a new compu- 
tational scheme which may be considered as an alternative to the ordinary delta expansion We introduce S as 
the parameter to dilate the scaling region of a given lattice model with the replacement of the lattice spacing a to 
(1 — <5) 1 / 2 a. As long as 5 can be tuned to some values close to unity, the lattice spacing a may be kept large enough 
in calculating physical quantities near the continuum limit. Further we perform the expansion in 5 to finite orders 
and setting S = 1 in the end of the calculation. Thus, our approach has some similarities with the ordinary delta 
expansion, and we use the term "delta expansion" to refer our method. We emphasize that our delta expansion on 
the lattice needs no extra parameter and the principle of minimum sensitivity [3] , both of which play, in the ordinary 
delta expansion, important but somewhat artificial roles to produce non-trivial results. 

To investigate and explore the above idea, we apply our method to the d = 2 non-linear a model on the lattice, 
which model is also called iV-vector model. In the present paper, we focus on two extreme cases, N = oo and N = 1 
which case corresponds to the Ising model. The model can be exactly solved in the large N limit, so we can examine 
to what extent our proposal is effective both in the qualitative and quantitative sense. The N = 1 case is also of 
interest, since the Ising model at d = 2 serves us a good testing ground of analyzing phase transition at non-zero 
temperature. 

II. DELTA EXPANSION IN SIMPLE EXAMPLES 

To illustrate our strategy, we first study two simple examples. 

A. Example 1 

Consider the problem of approximating the value of 0(0) where Q(x) is given as a finite series in I fx; 
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We introduce S by the replacement of x to x(l — 5) (0 < 6 < 1) in ^(a;), resulting a new function of two variables S 
and x, f2(x(l — 8)). Note that we can make the region around x = of f2(a:) wider by setting the value of S as close 
to 1. Then, if we can construct a truncated series of £l(x(l — 5)) in 1/x and 8 and it is effective at x > x* ~ O(l), the 
series would provide us the information on O(0) as long as 8 ~ 1. Here, we stress that for the plan to work the factor 
1—8 must be expanded in 8. As a systematic expansion of f2(x(l — 6)) effective at large x and small S, we employ 
the ordinary expansion in both variables around 1/x = 8 = 0. Then n th order approximant of fl(x(l — <5)), denoted 
as £l n (x,S), is written formally by 



& n (x, S) = c + 



E 



where 2j = 1/a; and zi = (5. For example, when n = 3, we have 
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In general, it is efficient to obtain the above expansion in the following manner: Consider 
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n„(x) = - - -5 + • 

x a; u; 

First shift x to a;(l — <5) and expand Q n (x(l — 6)) in d to the relevant order. For example the term x~ r should be 



expanded to the order 8 
dilation, we then have 



r , giving x r 
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Sctting S — 1 which means the infinite 
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From now on we use the symbol D as the operation of delta expansion to the relevant order with setting 8 = 1 
understood. To the order n we therefore obtain 



n„(a;,l) = D 
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It is easy to see that ( 1 — —j — > (n — > oo) for x > 1/2. Thus, when x > 1/2, 

n„(x,i) -» i = n(o). 

FIG. 1 shows the function 1) at n = 2, 3, • • • , 10. At those orders, there exists a plateu and the region of the 

plateu grows wider as the order increases. The value at a flat point which is a typical value on the plateu agrees with 
the value, S1(0) = 1. Thus the result is quite satisfactory. 



n=10 




FIG. 1: Plots of Q n (x, 1) at n = 2, 3, • ■ • , 10. The dotted line represents the value fi(0) = 1. 
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B. Example 2 

Next example deals with the relation of observables appearing in the d = 1 Ising model. In the model, the inverse 
temperature (5 is related to A/, the square of the screening mass (in lattice units) in the momentum representation, 
by 



1 / 4 
^=4 l0 H 1+ M 



At very low temperature the mass M is very small and we have the logarithmic relation 
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while at high temperature, M»l and 
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(2.3) 



Though the previous example has the limit, lim^^o £2(x) = 1, in the present example, j3 diverges logarithmically in 
the M — > limit. We show, however, that the delta expansion on the "high temperature expansion" (12. 3|) recovers 
numerically the asymptotic behavior of /3 near M ~ represented by l|2.2|) . 




FIG. 2: Behaviors of the "high temperature series" of (3 from 3 rd to 10 th orders. The dotted line represents the leading small 
M behavior. 



The series (|2 . 3[) itself breaks down at M = 4 (see FIG. 2) and it cannot be used for studying the small M behavior 
of (3. To improve the status, we implement the dilation of small M region around M = by shifting M — * (1 — S)M 
in (3(M). Using (|2.1j> . we then perform delta expansion at large M to give 
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To examine whether the above series captures the scaling behavior, we study the modification of the small M behavior 
due to the delta expansion. 

At small M, /3 is expanded as j3 = \ log jj + — + •■■)■ The leading term changes as 
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Truncating at 5 n and setting 5 = 1, we have 
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Note that as n — * oo, the constant part Y^k=i I diverges as logn. This reflects the logarithmic divergence of (3(M) 
in the M — > limit. The corrections by terms of positive integer powers experience drastic change: When the delta 
expansion is performed to a large order, lower order terms vanish when the value of 5 is tuned to 1. For example, we 
find that m — > m(l —S), m 2 — > rn 2 (l — 2<5 + <5 2 ) and m 3 — > m 3 (l — 3(5 H- 3<5 2 — <5 3 ), and these terms vanish themselves 



4 



at orders large enough. At small M, however, the delta expansion has a subtlety on the definition of the full order 
that is in accordance with the definition at large M . Though we do not know a convincing definition, let us proceed 
by supposing that we should confine ourselves with only the leading term at small M and it should be expanded in 5 
to the same order with the full order at large M. Thus, at n th order of delta expansion, ft n behaves at small M as 



\ fc=i / 



(2.4) 



As is obvious from the plot of ft n (see FIG. 3), the logarithmic behavior (|2.4p is recovered by the S- expanded 1/M 
series already at n = 3. Note that the scaling region develops to larger M as the order of expansion grows. In the 
scaling region to be seen in FIG. 3, small M behavior of ft{M) represented by (|2.2p is well approximated by subtracting 
j r from ft. Thus, even the limit of sequence {ft n } does not exist for any M as n — > oo, we can reproduce the 
scaling behavior both in the qualitative and quantitative respects. 
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FIG. 3: Behaviors of the 5-expanded 1/M series of ft at 3 rd and 10 th orders. The dotted lines represent the leading small M 
behaviors at respective orders. 



III. THE NON-LINEAR a MODEL AT LARGE N 



A. Brief review of the model 



The non-linear a model at two-dimensional Euclidean space is defined by the action, 



where the fields <J l (x) (i = 1, 2, • • • , N) obey the constraint, 

N 



a 2 {x) = ^2a i (x)a i {x)=N. 



The discretized space we work with is the periodic square lattice of the lattice spacing a where sites are numbered by 
two integers, (rii,n 2 ) = n. On the lattice the action may be written as 



s = 2ftj2°l- pYI ° n ■ ?n+e " - 

n n \i— 1,2 

where ft is defined as the inverse of the bare coupling constant /, 

1 



(3.1) 







and a n + e stands for the nearest neighbour spin of a n with ei = (1,0) and e 2 = (0, 1). The constraint is same as 
that in the continuum case, (?„ = N, and the first term in (|3.1[) is actually a constant that can be omitted. 
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Consider the correlation of fields, 

f3(a ■ <r r 



N 



7r / a 



d 2 p 



exp(— ip ■ ra)G(p, (3, a). 



By the use of Fourier transform method, one can calculate G(p, /3, a) to higher orders in j3. In particular, two moments 
S n ( cr oO'n) and J2 a ^W^n) are calculated to 21 th order [f|. One can use the result to obtain the correlation length 
£ as a series in j3. In the large N limit, however, it is more convenient for us to utilize the constraint and the fact 
that the structure of G becomes simple due to the absence of the wave function renormalization. Namely, taking the 
constraint a 2 = N into account and setting r = in the correlation (5q ■ a r ), we have 



d 2 p 
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U (27r) 2 M + 2E M= i )2 (l-cos^) 
Here we have rescaled the momentum by p — ► p/a and defined M by 

M 



(3.2) 
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where m represents the physical mass. Note that M is related to £ by M = £~ 2 . Since we consider in the next 
subsection the dilation of the region of M around M = 0, we express the inverse coupling (3 as a function of M. 
When the lattice spacing is small enough where M <C 1, we find from (|3.2|) 



Keeping only the first term, we have 
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(3.3) 



to 2 = 32a" 2 exp( — — ) = 32A|. 

This represents the dynamical mass in terms of the mass scale defined on the lattice. On the otherhand, when 
the lattice spacing is large where M 1, straightforward expansion of the right hand side of (|3.2p in M _1 gives the 
following: 



1 4 20 

M ~ U 2 + W 

From FIG. 4, it is apparent that the series 



112 676 
M 1 + IP 



4304 28496 



M 6 M 7 
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0(M~ S ). (3.4) 
~ 2 and the scaling behavior, (3 ~ 




FIG. 4: The graphs of j3{= 1//) as the function of logM. The dashed line represents the asymptotic behavior at small enough 
M, /? = — log || . Solid curves represent the graphs in the large M expansion from 2 nd to 9 th orders. 



log || , is not observed as it would. 



B. Delta expansion 

We perform the delta expansion of f3(M) with respect to a 2 or M following the manner adopted in example 2 
studied in section II. We shall see the asymptotically free behavior in the large M expansion of f3 and estimate the 
value of the non-pcrturbativc dynamical mass m in units of the lattice lambda parameter. 



6 



From (12.11). we find 



1 M \ 2 ) AP + { 3 ) M 3 



const 
n I M n 
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At small M, the leading term is transformed to 



1 M 
— log — 

4tt 5 32 



4tt V 5 32 ^ A: 
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We examine whether the scaling behavior written above is seen in (|3.5[) or not. FIG. 5 shows the plots of (3 n at 
n = 4, 9, 15 and 20. Dashed lines represent the leading small M behaviors at respective orders. We see that j3 n {M) in 
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FIG. 5: Graphs of /3„ at 4 th , 9 th , 15 th and 20 th orders as the function of logM. The dashed lines represent the leading 
asymptotic behaviors at small M at respective orders. 



1/M expansion is effective to log M ~ 1. In addition, we find that the logarithmic scaling behavior, (3 n ~ — -L log M, 
is seen at several and higher orders around logM ~ 2. Due to the dilation, the scaling region realized in 1/M 
expansion (|3.5p may develop toward larger M region. The tendency is confirmed from FIG. 5. 

Having observed the asymptotic scaling behavior in 1/M expansion, we can estimate the mass m in terms of 
even when the knowledge on the detailes of the small M behavior such as the information of the constant log 32 is 
absent. Let us write the scaling behavior as 

/?~--L(k> g M + C). (3.6) 

47T 

Note that (|3.6[) is derived only from the ultraviolet structure of the model. The value of the constant C(= — log 32) 
is, however, not obtained by the renormalization group argument alone. Then we like to show that our approach 
enables one to obtain the approximate value of the constant. 

We look for the matching point where the asymptotic behavior is supposed to begin and from there extrapolate 
fl{M) (at large M) to small M region. The matching point may be fixed by requiring the agreement of the derivative 
of P at large M with that of the leading term at small M which comes from renormalization group argument; 

d0(M)\ largeM = d 0(M)\ acaling _ J_ 
d log M dlogM 4tt' 
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Up to 20 th orders, the solutions, 
value of C can be estimated by 



M 



I, 4.7539, 4.6479, 



4.3597 exist at n = 4, 6, 8, ■ • • , 20 respectively. The 
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where M* denotes the solution. The result is —3.3066, 
the extrapolated asymptotic behavior, for example at n 



-3.3643,-3.3918, 
: 4, predicts 



-3.4383 at n = 4, 6, 8, • • • , 20. Then, 







scaling 



- — (log M- log 27.29). 
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(3.7) 



Comparison of (|3.7p to (|3.6p with C = — log 32 tells us that we have log 2 ^ 9 for log || . Now the approximation of 
the logarithmic constant leads us to the approximation of the dynamical mass. The result is m 2 as 27.29A^ at 4 
order. In the same manner, we obtain the ratio m 2 / h? L = 28.91, 29.72, 30.20, • • • , 31.13 at n = 6, 8, 10, 
values are obtained at higher orders and the sequence implies the convergence to the exact value. 
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C. Symanzik improvement 

In this subsection, we consider how to accerelate the speed of approaching to the asymptotic scaling and raise the 
accuracy of estimating the dynamical mass. 

The clue to the resolution is to note the presence of the subleading logarithmic terms, M k logM, in (3(M) at small 
M (see (|3.3I0 . The second and higher order contributions in (|3 . 3|) may prevent from the dominance of the leading 
term as we can see below: Consider the effect of these subleading logs in the application of the delta expansion ton 
P(M) at small M. If we expand powers and logarithms of M to S n , we find 

D[l] = 1, D[M r ] =0 (1 < r < n) 

and 

" 1 M M 2 
D[logM] =logM-Vr= D[M\ogM] = , D[M 2 \ogM] = 2! — — , 

nj ft Ti 

fc=i 

It is important to note that, though the terms of positive integer powers vanish, the logarithmic corrections survive 
after the delta expansion. Then suppose that M is large enough so that we can neglect the problem that to which 
order the first few terms should be expanded in 5. We then have 



k=l 



The third term of order M delays the scaling of (3 at finite n. 

The origin of the subleading logarithmic corrections is the propagator modified on the lattice. Actually, the 
expansion of the propagator at small p 2 reads 



1 1 1 I> 4 

M + T^Pl-TiT^P i + --- M + E^p* 12(M + £^ 



2 ^2 



(3.8) 



and the momentum integration yields —j^ log M from the first term and M log M from the second term. In general, 
k th term gives const. M fc_1 logM + regular terms in M. Thus, it is highly expected that Symanzik improvement @ 
accerelates the quick dominance since it subtracts the higher order p 2 corrections in the propagator. 

At the first order of Symanzik imporovement scheme, the coupling of spins at next-to-the nearest neighbour sites 
in both directions are introduced. It is well known that the resulting action becomes 

s = P 5^" ~ 3 XI ° n ' ° n+e x + Y2 X ^ n ' ° n+2e " 

n L |/ fi 
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Also in the improved action, we have the following result of the constraint at large TV, 

n7r d 2 p 1 



ft 



K (2tt) 2 M + 5 - § E M cosp M + i E M cos 2p M ' 



(3.9) 



Expansion of the denominator of the propagator for small p^ gives 

which has no p 4 contribution and thus MlogM term is absent in ft at small M. From (|3.9p . we find that ft behaves 
at small M as 

ft = -~t-(\og AI + C) + O(M), 

where C = —2.87298. This should be compared with (|3.6p . Though logM is invariant, the constant part which is 
not of universal nature is modified from C to C = C + 0.592754. The next-to-the leading logarithmic term, M log M, 
is absent and the approach to the asymptotic scaling would become faster than before because the correction to the 
scaling of ft is at most 0(M 2 ). 

From (|3.9p we find the large M expansion of ft, 



ft 
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1157 8419 3190499 
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M M 2 36M 3 36M 4 1728Af 5 
As is clear from FIG. 6, the above series is valid only for large M as in the case of the previous series (|3.4[) . However, 




FIG. 6: The graphs of the large M series of /3 in the first order Symanzik's improved action from 2 nd to 9 th orders 

once delta expansion is applied, we find that the large M series exhibits the correct logarithmic scaling behavior 
also for the improved action (see FIG. 7). Moreover, the extrapolated scaling function produces the following good 
approximate values for C = —2.87298, 



-2.8359, 



-2.8555, 



-2.8629, 



-2.8713, 



for orders n = 4, 6, 8, • • • , 20, respectively. The accuracy is much improved from those of the original action. 

To study further the results of Symanzik's improvement, we proceed to the second order. By introducing the 
spin-spin coupling of the form cr n • iJn+Se^ with the suitable weight, we can eliminate the second logarithmic correction 
M 2 log M in small M expansion of ft. The constraint equation becomes 



ft 



v d 2 p 



1 



\^ (2tt) 2 M + f - | £ M cosp M + ^ J2„ cos 2 Pfl - ^ E M cos 3 Pfi 
and the scaling behavior of ft is obtained as 

ft = -^(log M + C") + 0(M), 



where C" 



-2.72121 = C + 0.744526. 
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FIG. 7: Graphs of f3 n at 4 th , 9 th , 15 th and 20 th orders as a function of log M in the first order Symanzik action. The dashed 
lines represent the leading asymptotic behaviors at small M at respective orders. 



Now the large M expansion of /3(M) gives 

I 49 313733 2288857 727664156617 n(M -6 ] 
1 ~ M ~ 9M 2 + 8100M 3 ~ 7290A/ 4 + 262440000A/ 5 ~ ^ 

As in the case of the first order improvement, the scaling behavior is observed in the (^-expanded large M series and 
the extrapolation gives the approximation of the constant C" = —2.72121. At n = 4, 6, 8, • • • , 20, the results are as 
follows: 

-2.70458, -2.71498, -2.71823, -2.71956, -2.72097. 

The accuracy was further improved. We conclude that Symanzik's action plays a crucial role in the quantitative 
improvement on the delta expansion approach. 



IV. ISING MODEL AT d = 2 



Next we turn to discuss the Ising model at d — 2 which corresponds to N = 1 case of the iV-vector model. As 
well known, the second order phase transition at non-zero temperature is driven by the correlation length grown to 
infinitely large. Therefore the transition may be analyzed by dilation and delta expansion around the massless limit. 

From the calculation of the correlation function, the correlation length £ or M = £~ 2 would be obtained as a 
function of /?, the inverse temperature. In the Ising case, X^n^oCn) and ^ n n 2 (croCn) were calculated in Q to 25 th 
order in (3. We use the result and find that 

1 „ 5/3 67/3 3 38/3 5 9697/? 7 „„„ 8 1035164/? 9 352/3 10 

M = 2 + — + — — H — H 16/3 8 H 1 . 

2/3 3 45 189 4725 18711 3 

By inverting M~ l and j3, we obtain 

1 1 29 13 1503 1373 40581 9461 4553741 4309333 
^ ~ 2M ~ AP + 12M 3 ~ 2/\/ 4 + 80Af 5 ~ 24A/ 6 + 224M 7 ~ 16A/ 8 + 2304A/ 9 ~ 6407\/ 10 + " ' ' 



Now, near the transition point, conventional scaling form reads that 

h 



£ ~ const. I -rr — 1 
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where the critical exponent v and the inverse of the critical temperature (3 C are known to have values, v = 1 and 
P c = | log(l + V5) @|. Using £ 2 = M~ l , we then have 

P~ p c -AM&. (4.2) 

From (|4.2|) . we see that the present case has the limit limM^o P{M) = (3 C . However, the derivative of the first 
correction diverges as M — ► since the power of M, ^- is smaller than 1. This makes the convergence of P(M) to (3 C 
slow and forces us to have 1/M series of /3 to very large orders for the precise evaluation of f3 c (see FIG. 8). Then, 




02468 02468 



M M 

FIG. 8: (3 and fJ as functions of M at 25 th order. The dotted lines represent (3 C = \ log(l + sf2). 



rather than p c , we turn to the estimation of the critical exponent v . For the purpose, we consider g ^ M log[— g^ M ] 
of which behavior at scaling region reads 



4>{M) 







d log M 



log[ 



d/3 



d log M 1 2v' 
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The leading term of 4>{M) explicitly written in the right hand side of ()4.3|) comes from the second term of ([4.2 
Since the leading term of (f> is independent of M and invariant under the delta expansion, we expect that <f> - 
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d log M 



log[ 



d/3 



d log M 



] also behaves at small Mas^~^. Here P is given by 



Pn 



n\ 1 
1 ) 2M 



n \ 1 

2 M 2 



n\ 29 
3 J 12M 3 



n\ 13 
4 J 2M I: 
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FIG. 9 shows the plots of (j)(M) and 4>(M) in 1/M expansion at n = 25. The original function <f> plotted in the left 
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FIG. 9: Plots of <f>(M) and 4>(M) as functions of M at 25 th order. The dotted lines represent i(= \) 



side is far from the scaling region. On the other hand, ^-expanded series </> plotted in the right side shows greatly 
improved behavior. First we find that, apart from P(M) where the effects of the corrections to the leading constant 
P c is reduced but non-negligible, the corrections to ^ in <f> is surpressed enough by delta expansion and the stationary 
behavior is observed. We find that the scaling roughly starts about M ~ 10 and abruptly ends around M ~ 2.5 which 
signals the lower end of the region where <j) in 1/M expansion is valid. We can say that the asymptotic realm just 
starts at M ~ 5 and p n in 1/M expansion at n = 25 surely shows the expected behavior and strongly suggests the 
correct value of v. It is shown that the scaling region is enlarged as to be captured in 1/M expansion of p. 
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V. CONCLUSION 



Under dilation of the scaling region by shifting M — ► (1 — 5)M and the associated expansion in S, we have shown 
that 1/M expansion effective at large lattice spacing recovered the asymptotic scaling behavior in the non- linear a 
model at N ~ oo. In the large N case, Symanzik's improved action played an important role in the enhancement of 
the scaling. 

At N = 1, the scaling behavior was confirmed and rough estimation of the critical exponent v is possible. However, 
the precision is not so high even at 25 th order. Symanzik's improvement scheme would be an promissing candidate 
to raise the accuracy since, as shown in the N = oo case , the scheme would reduce the corrections to the asymptotic 
scaling. 

The present work stimulates us to investigate similar subjects on the lattice related to the approximation of the 
continuum field theories. Especially, it is of interest to apply the dilation with the delta expansion to lattice non- 
Abclian gauge theories. We hope to report the results in the near future. 
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